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ABSTRACT 

We present radial velocities for 85 cool stars projected onto the central parsec 
of the Galaxy. The majority of these velocities have relative errors of ~1 kms' 1 , 
or a factor of ~30-100 smaller than those previously obtained with proper mo- 
tion or other radial velocity measurements for a similar stellar sample. The error 
in a typical individual stellar velocity, including all sources of uncertainty, is 
1.7 kins -1 . Two similar data sets were obtained one month apart, and the total 
error in the relative velocities is 0.80 kms~ l in the case where an object is com- 
mon to both data sets. The data are used to characterize the velocity distribution 
of the old population in the Galctic Center. We find that the stars have a Gaus- 
sian velocity distribution with a mean heliocentric velocity of — 10.1±11.0 kms^ 1 
(blueshifted) and a standard deviation of 100.9±7.7 kms~ l \ the mean velocity 
of the sample is consistent with no bulk line-of-sight motion with respect to the 
Local Standard of Rest. At the 1 sigma level, the data are consistent with a 
symmetric velocity distribution about any arbitrary axis in the plane of the sky. 
We find evidence for a flattening in the distribution of late-type stars within a 
radius of ~0.4 pc, and infer a volume density distribution of r~ 1//4 in this region. 
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Finally, we establish a first epoch of radial velocity measurements which can be 
compared to subsequent epochs to measure small accelerations (1 kms^ 1 yr -1 ), 
corresponding to the magnitude expected over a timespan of several years for 
stars nearest to Sgr A*. 

Subject headings: Galaxy: center — techniques: spectroscopic — infrared: stars 



1. Introduction 

Evidence supporting the presence of a supermassive black hole in our Galactic Center 
(GC) is now very strong (Ghez et al. 2000; Schodel et al. 2002), and is certainly the strongest 
for any such object. While the current estimates of the central dark mass are nearly identical 
to the original estimates, M 3(1O 6 )M (Lacy et al. 1982; Serabyn & Lacy 1985), estimates 
of the bounding volume containing this mass have grown much smaller over the years. The 
resultant density estimates have grown from 10 6 M Q pc~ 3 to >10 17 M Q pc~ 3 , suggesting that 
the central dark mass is indeed a black hole. All of the mass determinations rely upon gas 
or stellar velocity measurements, either along the line of sight or in the plane of the sky, and 
are valid in the case that the motions are dominated by gravitational forces. 

The ionized gas in the central parsec orbits the black hole, and its bulk motion can be 
described by a model that assumes a large central mass (Serabyn & Lacy 1985; Giisten et 
al. 1987; Serabyn et al. 1988; Lacy et al. 1991; Roberts k Goss 1993; Herbst et al. 1993). 
The stellar motions in the central parsec exhibit ordered rotation for the young population 
(Genzel et al. 2000; Paumard et al. 2001; Genzel et al. 2003) and an isotropic distribution 
for the old population (Genzel et al. 2000). McGinn, Sellgren, Becklin, & Hall (1989) and 
Sellgren et al. (1990) use the line-of-sight motions of the old late-type stars to infer presence 
of the central dark mass; simliar experiments were done by Haller et al. (1996a) and Saha, 
Bicknell, & McGregor (1996). Finally, two groups have independently developed proper 
motion analyses that provide the best evidence supporting a super-massive black hole in the 
GC (Genzel et al. 1996, 1997; Ghez et al. 1998; Genzel et al. 2000; Ghez et al. 2000; Eckart, 
Genzel, Ott, & Sch6del 2002; Schddel et al. 2002; Ghez et al. 2003; Genzel et al. 2003). 



1 Data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific 
partnership among the California Institute of Technology, the University of California and the National 
Aeronautics and Space Administration. The Observatory was made possible by the generous financial 
support of the W.M. Keck Foundation. 
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Supermassive black holes have been inferred to exist at the dynamical centers of a 
number of galaxies having bulges, or bulge-like populations, and their presence may be ubiq- 
uitous in such galaxies (Ferrarese & Merritt 2000; Tremaine et al. 2002). It is important to 
establish the existence, and measure the properties, of supermassive black holes, given their 
potentially prominent role in controlling the energy output of galaxies and their diagnostic 
value in probing galactic evolution. In addition, it is also important to determine whether 
objects other than a supermassive black hole might contribute to enclosed mass estimates, 
i.e. various forms of dark matter. 

In this paper, we report the first data set obtained from a new observational program 
to determine the mass distribution in the central parsec, as determined from line-of-sight 
velocities of late-type stars. In this program, we obtained high spectral resolution data of 85 
cool stars throughout the central parsec using NIRSPEC on Keck II. We intend to combine 
these and similar data obtained in the future in order to measure accelerations, thereby 
constraining the mass distribution in the central parsec. 

2. Observations and Data Reduction 
2.1. Observations 

The observations were obtained on June 4, 1999, and July 4, 1999 with NIRSPEC, the 
facility near- infrared spectrometer, on the Keck II telescope (McLean et al. 1998), in high 
resolution mode covering i^-band wavelengths (1.98 /im to 2.32 /im). The long slit (24") was 
used in a north-south orientation, and the telescope was offset by a fraction of a slit width 
to the west between exposures. Two datacubes were built from the spectra obtained on the 
two nights. The field observed in July is almost entirely contained within the field imaged in 
June. Spectra of 40 of the same stars, 22 of them cool, appear in both data cubes. A log of 
observations is given in Table 1. The slit viewing camera (SCAM) was used to obtain images 
simultaneously with the spectra. These images make it easy to determine the slit orientation 
on the sky when the spectra were obtained; Figures la and lb show the inferred slit positions 
during all spectral exposures. From the SCAM images, we estimate seeing (FWHM) of / /6 
on the first night and Of! A on the second night. The plate scales for spectrometer and imager 
(0 / .'18) were taken from Figer et al. (2000). We chose to use the 5-pixel-wide slit (0'.'72) for 
the first slit scan and 3-pixel-wide slit (C'43) for the second slit scan in order to match the 
seeing on the respective nights. The corresponding resolving power was R~14,000 in June, 
and R~23,300 in July, measured from sky OH and arc lamp lines of known wavelengths in 
the two datasets. The slit positions were nearly parallel in the June scan, aligned with the 
long axis along the north-south direction, while they were quite skewed in the latter part 
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of the July scan because of problems with telescope tracking and/or image rotator control 
during these commissioning runs. 

The NIRSPEC cross-disperser and the NIRSPEC-7 and NIRSPEC-6 filters were used 
to image six echelle orders onto the 1024 2 -pixel InSb detector field of view. The approximate 
spectral ranges covered in these orders are listed in Table 2. Note that the two filters were 
used in concert on the second night in order to ensure complete rejection of light in unwanted 
cross-disperser orders. Some prominent features in our data set include diffuse and stellar 
lines of: Br-7 emission at 2.1661 /im in order 35, Fe III at 2.2178 /iin in order 34 of the June 
data, He I at 2.0581 jum in order 37, and the edge of the CO(2-l) transition at 2.2935 /im in 
order 33. 

Quintuplet Star ^3 (hereafter "Q3"), which is featureless in this spectral region (Figer 
et al. 1998), was observed as a telluric standard (Moneti et al. 1994). Arc lamps containing 
Ar, Ne, Kr, and Xe, and spectra of the sky, were observed to set the wavelength scale. In 
addition, a continuum lamp was observed through an etalon filter in order to produce an 
accurate wavelength scale inbetween arc lamp and sky lines (predominantly from OH). A 
field relatively devoid of stars ("dark spot") (RA 17 h 44 m 49!8, DEC -28° 54 6"8 , J2000) 
was observed to provide a dark current plus bias plus background image. A quartz-tungsten- 
halogen (QTH) lamp was observed to provide a "flat" image. 

2.2. Data Reduction 

We removed sky emission, dark current, and residual bias, by subtracting images of the 
dark spot. Due to the changing sky levels during the scan, the sky in the dark spot was 
scaled before subtracting it from target images. The scale factor was chosen such that it 
produced the minimum variance of residual sky line features in selected regions of the image 
after flat-fielding. The images were then "flattened" by dividing them by the QTH lamp 
image. These flattened images were then cleaned of bad pixels using a two-pass procedure. 
First, we replaced any pixel with a value 5a above neighbors within a 5 by 5 pixel area by 
the median of its neighbors in that area. Second, we replaced pixels with values that were 
higher or lower than the values of both immediate neighbors in the dispersion direction, 
and deviated by more than 10 times the poisson noise from those two neighbors. The 
results of these procedures are shown in Figure 2. The image was obtained with the slit 
positioned approximately north-south and centered on SgrA*. The echelle order number 
increases toward the bottom of the image, and north on the sky is downward in the slit. 
One can clearly see the diffuse Br-7 near the center of the image, and even the ionized "tail" 
associated with IRS7 (Yusef-Zadeh & Morris 1991). 
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Several images of the telluric standard, Q3, were taken in the same setup that was used 
to obtain the target images. In each of the images, spectra were imaged onto different rows of 
the detector. We then removed dark current, sky emission, and residual bias, by subtracting 
one image from the other. These images were further reduced by flat-fielding and cleaning 
deviant pixels, as described above for the target images. 

The etalon, arc and sky images were also reduced as described above. The reduced 
versions of these images were used for preliminary rectification of the warped echelle orders. 

2.3. Source Identification 

Figure 3 shows a SCAM image of the slit-scan region overplotted with identifiers taken 
from Table 4. Sources in the sample were identified based on their distance from the closest, 
brightest source in Genzel et al. (2000). Note that the sample in Genzel et al. (2000) is 
drawn from the sample in Genzel et al. (1996). In cases where no obvious counterpart in 
Genzel et al. (2000) could be found, we label the star according to the ID number in Table 4. 
All 85 cool stars in the table are labelled in the figure, and the spectroscopic sample is ~50% 
complete to K=13. 

2.4. Spectral Image Rectification 

We rectified the data format into orthogonal axes of wavelength and (approximately) 
declination by mapping a set of observed features in the warped images onto a set of dewarped 
grid points. Typically, 15 to 20 spectral lines of known wavelength were used in each order 
as wavelength fiducials; the zero-point wavelength was set by sky OH lines in the target 
spectra. Two stellar continuum spectra and two slit edges traced across the length of the 
dispersion direction served as spatial references. Figure 4 shows the result of this process, 
which we explain in more detail below. 

In general, we traced arc, sky, and etalon lines as follows. Each line was divided into 
10 to 20 equally sized samples along the slit length, each sample containing 10 rows of data. 
Rows within each sample were averaged together, and the location of the centroid along the 
x-axis was recorded. We then used a 3rd order polynomial to fit the location of the centroids 
as a function of location along the y axis. In cases where the centroiding was less robust, due 
to a nearby unidentified spectral line, or when remaining bad pixels affected the centroiding, 
the order of the polynomial fit was reduced and/or points near the edges were deleted. We 
traced spatial locations along the slit length as a function of location along the wavelength 
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axis in much the same way, except that the centroids were computed along the y axis. 

Given the spatial and spectral centroids from the preceeding process, we then rectified 
the data in three stages: 1) rectify images using the arc and sky lines, 2) measure the 
etalon line wavelengths in the rectified image, and obtain a solution to the etalon equation, 
and 3) use the analytically-determined wavelengths of the etalon lines to produce a better 
rectification matrix. 

In the first stage, we produced a solution using an order containing many (15 to 20) arc 
and sky lines. Order 35 of the June dataset and order 36 of the July dataset were selected. 
OH line wavelengths were taken from Abrams et al. (1994). The arc line wavelengths were 
obtained from the National Institute of Standards and Technology (NIST) Atomic Spectra 
Database 9 . The arrays of coefficients from spatial and spectral fits were then used to produce 
a mapping between points in the warped frame, and points in the dewarped frame. We then 
fit a two-dimensional polynomial, of second or third order, to these points, resulting in a list 
of transformation coefficients. Finally, we rectified the selected orders of the etalon image. 

In the second stage, we measured the wavelengths of the etalon lines using "SPLOT" 
in IRAF 10 . This provided preliminary wavelengths for each etalon line. Exact etalon wave- 
lengths are given by the etalon equation, i.e. A=t/2n, where t is the thickness of the air gap 
of the etalon and n is the etalon order number. Each etalon line will have its own equation 
that relates its specific wavelength to its specific order number. We simultaneously solved a 
set of etalon equations for the etalon lines we fit using the wavelength solution from stage 
one, under the constraint that all equations must use the same thickness, and that the order 
numbers be increasing integers from longer to shorter wavelengths. We varied the thickness 
and order numbers until the wavelengths predicted by our fit were closest to the wavelengths 
found in stage one. 

In stage three, we used our etalon wavelengths from stage two to fit the etalon lines 
in all echelle orders. Together with the same spatial information used to produce the first 
stage rectification matrix, a new rectification matrix was produced. The improved quality 
of the etalon lines allowed for a higher order (fourth or fifth order) polynomial fit to the 
etalon mapping between warped and dewarped points. The new matricies were applied to 
the appropriate spectral orders of target images to rectify them. 



9 http://physics. nist.gov/PhysRcfData/ASDl/nist-atomic-spectra. html 

10 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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2.5. Building the Data Cubes 

We extracted stellar spectra, typically four to eight per slit pointing, using "APALL" 
in IRAF. Spectra of the same object observed in multiple telescope pointings were coadded, 
after first shifting each spectrum to account for slight variations in the object's position 
along the slit width (dispersion) direction. This was done by first translating each spectrum 
along the wavelength direction by an amount equal to the offset between it and the spectrum 
of IRS7. To obtain this offset, we cross-correlated the spectra of the target and IRS7 near 
2.05 //in, where many telluric features are present. We applied another shift to account 
for the fact that IRS7 was not likely in the center of the slit when it was observed. To 
estimate this shift, we measured the median of the offsets of all target stars with respect 
to the spectrum of IRS7, and assumed that the distribution of positions within the slit was 
random, i.e. we added the negative of this offset to each target spectrum. The total shift to 
compensate for target position within the slit along the dispersion direction was typically a 
few tenths of a slit width, or equivalently, approximately one pixel. 

In using the procedure above, we are relying on the fact that wavelengths of the telluric 
features do not change. In addition, we rely on the fact that the grating remained in the same 
position throughout the observing sequence, and that variable flexure in the instrument was 
inconsequential; note that NIRSPEC sits on the Nasmyth deck and therefore has a constant 
orientation with respect to gravity. Our data do not suffer from the problem that currently 
plagues some NIRSPEC echelle data, that of a wavelength shift induced by grating motion 
when the image rotator is slewed. 

Spectra of the source Q3 were extracted and shifted to the center of the slit, median- 
combined, and then normalized to 1 in order 36. Target spectra were then divided by these 
final standard spectra to remove telluric lines and to normalize and remove the shape of 
the stellar continuum, making it easier to compare the CO bandhead and other features of 
interest. 

On the first night of observations, we obtained spectra for 152 stars, including 102 
cool stars, 34 hot stars, and 16 stars that have featureless or indeterminate spectra. We 
have separated the cool stars into two groups, based upon quality of the data, i.e. 78 are 
of high quality and 24 are of low quality. On the second night of observations, we obtained 
spectra for 68 stars, including 38 cool stars, 16 hot stars, and 14 stars that have featureless 
or indeterminate spectra. Of the cool stars, 29 have high quality spectra, and 9 have low 
quality spectra. 

In this paper, we present the velocities of only the cool stars, given that their spectra 
can be analyzed to produce the most precise velocities; in addition, the cool star orbits are 
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much more likely to be thermalized and thus better suited for use in probing the local matter 
distribution than the hot star orbits. We only include the stars with high quality spectra. 
Unfortunately, the hot star spectra usually have intrinsic features, produced in photospheres 
or winds, that are confused with features from nearby ionized gas, i.e. gas in the GC mini- 
spiral. In addition, their features are sometimes broad and/or asymmetric, both of which 
lead to much larger errors in determining their velocities compared to those obtainable from 
cool star spectra. We will report on our efforts to extract precise velocities from the hot star 
spectra in a separate paper. 



2.6. Cool star velocities 

Most (~60%) stars in our sample are cool, showing varying amounts of absorption 
due to CO, as well as many other atomic and molecular species. We used the CO bandhead 
(edge rest wavelength is 2.2935 /im) to determine velocities by cross-correlating spectra of the 
sample stars with vacuum rest frame spectrum of Arcturus (Hinkle et al. 1995) (R=100,000), 
correcting for Earth motion toward the Galactic Center and differences in spectral type 
between the sample stars and Arcturus. 

We corrected Arcturus' spectrum to the rest frame by maximizing the cross-correlated 
power between it and a spectrum of the Sun in the wavenumber region from 4300 cm' 1 to 
4375 cm -1 , after first resampling the solar spectrum (Livingston & Wallace 1991) so that 
both spectra had the same effective resolution (R~100,000). We believe that the adjustment 
to the rest frame is accurate to within one tenth of this resolution, or ±0.30 kms" 1 . The 
rest spectrum of Arcturus was then resampled and smoothed to match the spacing and 
resolution of the GC target spectrum. We measured geocentric velocities of the target stars 
by cross-correlating target and template spectra in the 2.29 fim to 2.31 /xm region. Figure 5 
shows a full spectrum of IRS7 for observed orders, and Figures 6 through 9 show sample 
spectra of target stars compared to the rest spectrum of Arcturus in the region covering the 
CO bandhead. 11 The spectra have been shifted along the wavelength scale to maximize the 
cross-correlated power between the target and template spectra. 

We adjusted the geocentric velocities to the heliocentric frame by removing Earth mo- 
tion, +6.6 kms~ l in June and —8.1 kms^ 1 in July, where a positive value indicates motion 
toward the GC. The error in the resultant velocity due to this adjustment is estimated to be 
±0.10 kms- 1 . 



11 A complete set of such figures for all stars in this paper can be found at: 
http://www.stsci.edu/ figer/papers.html. 
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A correction to the velocities was also made to account for the differences in the spectral 
types of target stars versus that of Arcturus (K1.5III). This correction was determined in a 
two-step process. First, we cross-correlated template spectra of giants in the Wallace-Hinkle 
(WH) atlas (Wallace & Hinkle 1996) with that of Arcturus, and assigned a velocity correction 
as a function of temperature, as shown in Figure 11a. We removed data for variable stars 
(open circles in the plot) because their measured velocities can be influenced by gas motion in 
their atmospheres. Second, we determined velocity corrections for the target stars based on 
the depths of the CO feautures in the spectra; the relation between the CO equivalent width, 
EWco, and temperature is shown in Figure lib. The velocity offsets and temperatures 
for the template stars were taken from Ramirez et al. (1998). For the equivalent width 
measurements, we defined the continuum as the average flux from 2.288 /im-2.293 /iin, and 
we isolated the region from 2.294 /xm-2.300 /zm for measuring EWco- Note that this process 
produces a velocity correction of nearly zero for a K1.5III star, i.e. Arcturus, as expected. 
The error in the velocity correction has contributions from the error in the measurement 
of EWcoj the error in the relation between temperature and EWco, and the error in the 
relation between the velocity correction and temperature. We assign an error of 2 A to the 
EWco, or 166 K in temperature. The error in the temperature versus EW C o relation is 81 K. 
Together, these two errors contribute 0.8 kms~ l error in the velocity correction. Summing 
this, in quadrature, with the error in the velocity correction versus temperature relation 
gives a total error of 1.39 kms' 1 . There is another potential source in error if the target 
stars are long-period variables. In that case, the outer atmospheres exhibit variable radial 
motions with respect to the systemic velocities of the stars. This might be an important 
consideration for the latest types in the target sample. 

In summary, we estimate the following systematic errors in determining heliocentric 
stellar velocities: Arcturus rest frame adjustment (±0.30 kms^ 1 ), Earth motion adjustment 
(±0.10 kms~ r ), and spectral type correction (±1.39 kms^ 1 ). The total systematic error is 
then ±1.43 kms" 1 ; note that the relative velocity errors are smaller and can be determined 
from the data (see below). 

Table 4 lists identifiers, relative positions, estimated spectral types, heliocentric veloci- 
ties, and internal velocity errors for each cool star, where V* = V mcasuTC d + V ffi — V corr . 

3. Analysis and Results 

In this section, we present an anlysis of the statistics concerning the velocity distribution 
and associated errors. 
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3.1. Statistics and Errors 

Table 5 gives statistics of the measurements presented in Table 4. In particular, we 
measured velocities for 78 stars in the first data set, and 29 stars in the second data set. 
There are 22 stars in common in both data sets, and there are 85 unique stars in the combined 
data set. Figure 12 shows the histogram of velocities from the combined data set, along with 
a similar plot for data from Genzel et al. (2000). Figure 13 shows a "bubble plot" of stars and 
their velocities from the combined data set. The unweighted average (— 10.1±11.0 kms^ 1 ) 
and median (— 7.5±11.0 kms~ l ) velocities for the combined data set are consistent with 
each other and consistent with the expected value resulting from the projection of the Sun's 
peculiar velocity with respect to the Local Standard of Rest (LSR), —10 ± 0.36 ferns -1 
(Dehnen & Binney 1998); this agreement suggests that there are no large systematic errors 
in the absolute velocity calibration. Strictly speaking, however, it is possible that the target 
population could exhibit bulk motion along the line of sight with respect to the LSR (Miller 
1996). 

The mean difference in velocities, in cases where we measured two velocities for a sin- 
gle star, is —0.80 ferns -1 , with a standard deviation of 1.53 ferns -1 . We assume that the 
individual errors are uncorrelated, so that the la velocity error for cases with two velocity 
measurements, is 0.77 ferns -1 , and 1.08 ferns -1 for cases with one velocity measurement. 
Thus, the errors in Table 4 are small, and the velocities are consistent with those in Genzel 
et al. (2000), as seen in Figures 12 and 14. The standard deviation of velocities for our com- 
bined data set is 100.9±7.7 ferns -1 , where we take the error to be dominated by statistics, 
i.e. error=cr/(2iV) a5 . 

3.2. Asymmetry 

We examined the data for asymmetric distributions in order to investigate the dynamical 
state of stars in the sample. To first order, the velocity distribution appears to be Gaussian 
(Figure 12). To inspect the velocity distribution further, we examined the data for rotation 
about axes along galactic longitude and latitude. Figure 15 shows the line-of-sight velocities 
(v) as a function of the Galactic longitude and latitude (the Sun-GC distance is assumed to 
be 8 kpc). It appears that there is no significant correlation between velocity and projected 
distance from Sgr A*, implying that the 85 sample stars do not show any considerable 
figure rotation. To explore this question in a more statistically robust way, we examined the 
observed velocity distribution for rotation about an arbitrary axis in the plane of the sky. 
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To do this, we defined an asymmetry statistic, 

A = [n rcd - ^bluc] aaxis<PA<Qax . s+180 o - ["-red - n bluc] aax . s+180< p A<aaxis+360 ° , (1) 

where "red" and "blue" refer to redshifted and blueshifted velocities, and the terms in the 
brackets refer to the number of stars with position angles between the angles given in the 
subscripts, where a ax i S is the position of an assumed axis. The position angles are measured 
East of North, for an axis with the origin at the location of Sgr A*. 

The normalized asymmetry statistic is shown in Figure 16a for all possible assumed 
position angles for the axis. The solid line has a maximum value at a position angle of 0°, 
a somewhat surprising, and perhaps coincidental, result. Otherwise, there appears to be no 
significant asymmetry by this statistic. We performed Monte Carlo simulations in order to 
determine the significance of the peak in A. 

Figure 16b shows the cumulative normalized probability of a model system having a 
normalized asymmetry statistic, A/N, greater than the value on the x-axis. Model systems 
have 85 stars, each with uniform random locations within a 1 by 1 pc box centered on Sgr A* 
and normal random velocities with <x=101 kms^ 1 . One thousand systems were generated, 
and A was determined for position angles between zero and 180° at 1.5° increments. The 
maximum A was then stored for each system. The curve shows that approximately 20% of 
systems will have A at least as great as 0.25, the value we observe. This suggests that the 
observed system is consistent with a symmetric distribution in the asymmetry statistic to 
within approximately la. 



4. Discussion 

In this section, we discuss the results and their implications for determining the dynam- 
ical properties of the late-type population in the GC. We also describe their correspondence 
with the data from the 198 cool stars in Genzel et al. (2000), 106 of which are located with 
the region covered by our observations, and 49 of which are common with stars in our data 
set. In addition, we discuss the possibility of using the data as a first epoch of observa- 
tions that will yield measurements of orbital accelerations upon combining data from future 
observations. 

Our results confirm that the sample population traces a distribution of stars in isotropic 
orbits around a large dark mass. While a few of the brightest stars in our sample are young, 
i.e. IRS7, they evidently do not influence the statistics enough to imprint on them any 
velocity asymmetry in the statistics we have examined. 



4.1. Comparison to Radial Velocities in Genzel et al. (2000) 



There are 49 stars in common between our data set and that in Genzel et al. (2000). For 
a detailed comparison between the two data sets, we have plotted the normalized asymmetry 
statistic and associated probability function in Figures 17a and 17b from data in Genzel et 
al. (2000). In general, the asymmetry curves in Figures 17a and 16a are similar, showing 
maxima near 0°, and two other broad humps at larger angles. Curiously, the Genzel et 
al. (2000) data show a conspicuous number of stars with positive velocities on either side 
of a rotation axis, regardless of position angle, as seen in the dotted and dashed lines. 
This is consistent with the fact that the median and mean velocity in their data set are 
+12.0 kms^ 1 and +2.2 kms^ 1 , with an error of 7.3 kms^ 1 , as opposed to —10.1 kms^ 1 
and —7.5 kms^ 1 , with an error of 11.0 kms^ 1 , in ours; both data sets are consistent with a 
Gaussian velocity distribution having similar standard deviations, 100.9+7.7 kms^ 1 in our 
data and 103+5.2 kms~ l in theirs. 



Both data sets also produce similar results when used to infer an enclosed mass, assuming 
a relationship between the velocity dispersion and the local gravitational potential, i.e. the 
Bahcall-Tremaine (Bahcall & Tremaine 1981) estimator: 



where, G is the gravitational constant, N is the number of stars, v is the line-of-sight velocity, 
R is the projected radius, w is the normalized statistical weight, and the sum is over all test 
particles within the enclosed projected radius, Ro. The individual weights are assigned to be 
one over the square of the errors in Table 4, i.e. the inverse variances, and the mean velocity 
of the sample is subtracted before applying Equation 2. 

Figures 18a, 18b, 19a, and 19b show the velocity dispersion, a, and enclosed mass, M B t, 
as calculated using data in Table 4 (a) and from Genzel et al. (2000) (b); we shifted velocities 
in both data sets to adjust the mean of the sample to zero. The error bars are the la values 
calculated from the data. The Genzel et al. (2000) data extend to smaller radii, but both 
data sets show the trend of a few million solar masses of enclosed material at the smallest 
radius, and a steady increase of enclosed mass for larger radii. Within errors, the data sets 
agree in their prediction of enclosed mass. Because it is tempting to over-interpret these 
figures, we limit our conclusions to be that the enclosed mass within a projected radius of 
0.1 pc is ~2±1(10 6 ) M according to the BT statistic. This estimate compares favorably 



4.2. 



Enclosed Mass Estimates 




(2) 



R<Ro 



13 



with others (Lacy et al. 1982; Serabyn & Lacy 1985; McGinn, Sellgren, Becklin, & Hall 
1989; Sellgren et al. 1990; Haller et al. 1996a; Genzel et al. 1997; Ghez et al. 1998; Genzel 
et al. 2000; Ghez et al. 2000; Eckart, Genzel, Ott, & Schodel 2002), although they tend to 
be low with respect to the "true" enclosed mass, 3.45±0.5(10 6 ) M & , inferred by assuming a 
Keplerian orbit for one of the closest stars to the dark mass (Schodel et al. 2002; Ghez et al. 
2003). 

One notable weakness of the BT estimator in this application is the fact that the dis- 
tribution of observed particles may deviate from that assumed in the derivation of the es- 
timator. In particular, our coverage represents a pencil-beam that slices through the GC, 
thereby including test points significantly farther from the Center than the projected extent 
of our sample in the plane of the sky; this will tend to add objects with lower line-of-sight 
velocities than would be the case if we were completely sampling all projeced radii in the 
system. In addition, our sample includes stars on plunging orbits, i.e. with high line-of-sight 
velocities, without containing the corresponding group of stars from the isotropic distribution 
at projected radii beyond what was sampled. This bias will tend to increase the estimated 
projected mass. Haller et al. (1996b) highlight these biases and conclude that one could not 
likely estimate the enclosed mass in the central parsec to better than a factor of two using 
the BT estimator. An additional bias can be introduced if particles are missing from the 
distribution, i.e. in the case of an evacuated cavity, as we shall see later. 



We attempted to fit the radial profiles and velocity dispersions assuming a relatively flat 
distribution, within some radius, r , that is described by a power law of index n, according 
to equation 3, and similar to the form used by Sellgren et al. (1990). 



We find good agreement between the models and observations for ro ~0.4 pc, and 
n~ —1/4, using the following procedure. We created a set of model test points in a Monte 
Carlo simulation assuming equation 3 from to 10 pc, and a total enclosed mass at the outer 
radius of 3(10 7 ) M Q (McGinn, Sellgren, Becklin, & Hall 1989). Each particle is governed 
by a gravitational potential set by the the enclosed mass within the radius of the particle, 
including the masses of particles with orbits at smaller radii and a central point mass to 
represent the black hole, having M B h=3.45(10 6 ) M , the mean of the values in Schoedel, 
Genzel, Ott, & Eckart (2003) and Ghez et al. (2003). Each particle is on a circular orbit, with 
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uniform random phase and inclination with respect to the line of sight. The total sample 
size ~10 6 particles. We used the simulation to find the most probable r of ~0.4 pc. We then 
determined the best-fit index by minimizing the of the difference between the model and 
observed velocity dispersion (Figures 21a and 21b). The results are shown in Figures 22a 
and 22b, where we see that the best fit is for a power law index ~ — 1/4 for r<r . 

Figures 23a and 23b show plots of locations for 10,000 randomly selected points in a 
simulation with n=— 2 (a) and one with n= — 1/4 (b). Figures 24a and 24b show the projected 
densities and velocity dispersions for the modelled particles, for circular and isotropic orbits. 
The surface density increases as 1/r outside or ro, as expected, and rolls over to a roughly 
constant value within this radius. The resultant velocity dispersion is relatively flat within 
r , consistent with the observed behavior. 

Sellgren et al. (1990) and Haller et al. (1996a) describe the absence of the CO bandhead 
feature in spectra of integrated light within 0.6 pc of the GC. They consider the hypothesis 
that CO in the atmospheres of late-type stars is destroyed by photo-dissociation or stellar 
collisions in the region. Others have also suggested yet other explanations for the CO hole, 
i.e. stellar mergers/collisions (Lacy et al. 1982; Lee et al. 1996; Alexander 1999; Bailey & 
Davies 1999), atmospheric stripping (Alexander & Livio 2001), and atmospheric heating 
(Alexander & Morris 2003). Genzel et al. (1996) provide strong evidence for a hole of late- 
type stars in the form of a steep decrease in the projected number density of such stars 
(their Figure 6) and a decrease in the strength of the CO bandhead feature in integrated 
light within 5" of the GC. They also invoke these features as evidence for the destruction of 
the largest red giants in interactions with main sequence stars. Our findings are consistent 
with previously obtained evidence for a flattening of the density distribution; however, we 
consider that yet another physical process may explain the observed phenomenon. Miralda- 
Escude & Gould (2000) argue that a cluster of ~24,000 stellar mass black holes should have 
migrated toward the central parsec through dynamical friction on a timescale of ~10 Gyr 
(Morris 1993). They find through simulations that this cluster will eject stars on a time 
scale of ~1 Gyr. We propose that such a cluster of black holes could provide the scattering 
sources for ejecting relatively low mass, and old, stars from the region. 

We favor this hypothesis versus the alternatives for the following reasons. First, the 
requisite star formation needed to produce the stellar black holes is supported by the presence 
of young and intermediate-age stars in the GC, i.e. Lebofsky, Rieke, & Tokunaga (1982), 
Rieke (1987), Lebofsky & Rieke (1987), Haller & Rieke (1989), Haller (1992), Rieke et 
al. (1993), Blum, DePoy, & Sellgren (1995), Krabbe et al. (1995), Genzel et al. (1996), 
Narayanan, Gould, & Depoy (1996), Sjouwerman et al. (1999), and Figer (2003). Second, 
Morris (1993) and Miralda-Escude & Gould (2000) predict that dynamical friction should 
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bring black hole remnants from massive stars into the central parsec in substantially less 
than 10 Gyr. Alternate hypotheses require interaction (mergers or collisions) rates that may 
be in excess of what is predicted by the observed stellar density distribution and/or stellar 
evolution scenarios that are poorly understood, i.e. the spectral energy distribution emitted 
by a post-merger star; however, see Genzel et al. (2003) for an alternate interpretation. 

Finally, we explore the difference in enclosed mass estimates from the proper motion 
studies, M~2.5-3.2(10 6 ) M (Genzel et al. 1996; Ghez et al. 1998), and those from more 
recent estimates derived from single orbits, M=3.45±1.5(10 6 ) M (Schoedel, Genzel, Ott, & 
Eckart 2003; Ghez et al. 2003). The former are inferred using a mass estimator, such as the 
BT estimator. We attribute the difference as due to differences in the observed distribution 
versus that assumed in the derivation of the BT estimator. To take an example, consider 
Figures 25a and 25b, in which we plot the true enclosed mass and that estimated by the BT 
estimator, versus radius, for a system with no hole (a), and with a hole having r =0.4 pc 
(b). In this case, we are using circular orbits, so the normalization factor in equation 2 is 
32/3, instead of 16. The plot shows that the BT estimator overestimates the enclosed mass 
at large radii, as expected from an incompletely populated system, and underestimates the 
enclosed mass at radii less than the hole radius. Note that the true mass and that inferred 
using the BT estimator converge for radii above r up to ~1 pc. This effect is also seen in 
the observations (Figure 19) where the enclosed mass estimates converge to ~3.8(10 6 ) M 
for a projected radius above 0.4 pc, in good agreement with the black hole mass as inferred 
from single orbits. 

We have used the models in this section to explore the plausibility that a hole in the 
distribution of old stars can reproduce the observed velocity dispersion and projected number 
distributions. While the models are consistent with the presence of such a hole, we look 
forward to more robust modeling that can be tested with statistical methods to determine 
the most likely volume density of old stars in the GC. If a hole is required, then it will be 
interesting to determine the physical mechanism that produces it. In particular, it will be 
interesting to determine whether the Chandra sources in the central parsec are consistent 
with the presence of a cluster of black holes (Muno et al. 2002). 

4.4. Accelerations 

We expect to be able to measure accelerations by comparing velocities in this paper to 
those collected in the future. The magnitudes of these accelerations depend on the projected 
radius and orbital inclination, as shown in Figure 26, where the potential is defined by a point 
mass having M=3.45(10 6 ) M Q . For an "average" case, we expect to measure an acceleration 
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that exceeds the errors in 5 years from the first epoch of observations in this paper. Note 
that the error in the estimated acceleration only depends on the statistical errors of the 
measurements, and not the systematic errors, i.e. the peculiar solar velocity with respect to 
the local standard of rest. 



5. Conclusions 

We have established a first epoch of radial velocity data for a large (85) number of 
cool stars in the central parsec of the Galaxy. The data set have much higher precision, 
by at least a factor of 30, than those previously obtained for similar samples. The data 
agree with previous results, and we conclude that the late-type population has Gaussian 
velocities in the line-of-sight component and that the enclosed mass at a projected radius of 
0.2 pc is consistent with a value ~3.45(10 6 ) M . We have identified dynamical evidence for 
a flattened distribution of late-type stars within the central parsec, with an edge at ~0.4 pc. 
The finding adds evidence of such a flattening, in addition to that from previous studies 
finding an absence of the CO bandhead feature in the integrated light of late-type stars in 
the region. While the integrated light measurements might suffer from dilution by the light 
of young stars, the dynamical evidence does not suffer from such a bias. 

We acknowledge useful conversations with Mike Fall. We also acknowledge the work 
of: Maryanne Angliongto, Oddvar Bendiksen, George Brims, Leah Buchholz, John Canfield, 
Kim Chin, Jonah Hare, Fred Lacayanga, Samuel B. Larson, Tim Liu, Nick Magnone, Gunnar 
Skulason, Michael Spencer, Jason Weiss and Woon Wong. In addition, we thank the CARA 
instrument specialist Thomas A. Bida, and all the CARA staff involved in the commissioning 
and integration of NIRSPEC. 
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Table 1. Log of Observations 



Name a 


Resolution 13 


Filter c 


Integ. 


Frames 


Slit Size 


Date 


GC1 


14,000 


NIRSPEC-7 


60 s. 


40 


0772 x 24" 


4 June 1999 


GC2 


23,300 


NIRSPEC-6/7 


150 s. 


67 


0743 x 24" 


4 July 1999 



a "GCl" refers to the first slit scan while "GC2" refers to the second slit scan. 
b The resolution is A/AAfwhm, where AAfwhm is the half-power line width of unre- 
solved arc lamp lines. The slit width was 5 pixels in GC1 and 3 pixels in GC2. 

c NIRSPEC-7 has half-power points of 1.85 /im and 2.62 //m. NIRSPEC-6 has half- 
power points of 1.56 /im and 2.30 pm. Half-power points are from Figer (1999c). 
Because the orders are longer than the width of the detector, the spectra are non- 
contiguous in wavelength. 

Note. - All images were obtained with the slit positioned approximately north- 
south. The multiple correlated read mode (Fowler & Gatley 1990) with 16 reads at 
the beginning and end of each integration was used for both data sets. 
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Table 2. Wavelength Coverage of Echelle Orders (jum) 



Echelle 


GCl a 


GCl a 


GC2 


GC2 


Order 


NIRSPEC-7 


NIRSPEC-7 


NIRSPEC-6/7 


NIRSPEC-6/7 




min 


max 


min 


max 


33 


2.281 


2.315 


2.290 


2.324 


34 


2.214 


2.248 


2.223 


2.256 


35 


2.152 


2.184 


2.160 


2.192 


36 


2.092 


2.124 


2.101 


2.131 


37 


2.036 


2.067 


2.044 


2.074 


38 


1.983 


2.013 


1.991 


2.020 



a "GCl" refers to the first slit scan while "GC2" refers to the second slit scan. 
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Table 3. Velocity Correction for Spectral Type 



Wri atlas 


Voff 


Type 


rp c 


V d 

v corr 


Star 


km s^ 1 




rv 


km S 


-l 


Kap Gem 


1.86 


G8III 


4866 


2 


31 


Xi Her 


3.93 


G8.5III 


4788 


1 


97 


lot Cep 


1.60 


K0-III 


4565 


1 


01 


11 Cep 


0.46 


K0.5III 


4495 





71 


Alp Boo 


-0.64 


K1.5IIIp 


4360 





13 


Kap Oph 


-0.46 


K2III 


4294 


-0 


15 


39 Cyg 


-2.00 


K2.5III 


4232 


-0 


42 


31 Lyn 


-1.89 


K4.5III 


4000 


-1 


42 


Gam Dra 


-2.25 


K5III 


3946 


-1 


65 


Alp Tau 


-1.06 


K5+III 


3946 


-1 


65 


Mu UMa 


-4.79 


M0III 


3845 


-2 


09 


Gam Sge 


-1.17 


M0-III 


3845 


-2 


09 


Nu Vir 


-1.94 


M1III 


3752 


—2 


49 


75 Cyg 


-1.65 


MHIIab 


3752 


-2 


49 


HR8989 


-1.61 


M2III 


3666 


-2 


86 


Chi Peg 


-2.42 


M2+III 


3666 


-2 


86 


HR8621 


-10.69 


M4+III 


3517 






RLyr 


-4.34 


M5III 


3454 






HR8530 


-14.91 


M6III 


3399 






EU Del 


-4.14 


M6III 


3399 






SW Vir 


-4.85 


M7III: 


3351 






BK Vir 


-10.58 


M7-III: 


3351 






RX Boo 


-4.47 


M7.5-8III 


3320 







a Giants from the Wallace-Hinkle atlas (Wallace 
& Hinkle 1996). 

b Measured velocity offset with respect to Arc- 
turus. 

c Temperatures from Ramirez et al. (1998). 
d Fit velocity offset with respect to Arcturus. 
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Table 4. Velocity Data 



Identifier 


Designation 11 


ARA b 


ADEC b 




Sp. c 


Vi d 


v 2 ° 


Vave' 




VGenzcl' 1 


"Gonzcl h 










" 


Type 


km s 1 


km s 1 


km s 1 


km s 1 


km s 1 


km s 1 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


259 




-0.68 


-1.72 


1.85 


K4III 


95.7 


97.4 


96.5 


0.8 






313 


IRS29S 


— 1.79 


0.82 


1.97 


K3III 


-160.3 




-160.3 


1.1 


-93 


20 


327 




-1.58 


1.24 


2.00 


G9III 


-156.0 




— 156.0 


1.1 






360 


G288? 


— 1.28 


2.61 


2.91 


G9III 




133.7 


133.7 


1.1 


17 


30 


365 




1.23 


2.79 


3.05 


K2III 


1.9 




1.9 


1.1 






221 




0.99 


-3.36 


3.50 


M6III 


-16.0 




-16.0 


1.1 






222 


IRS33W 


-0.03 


-3.51 


3.51 


K4III 


72.0 


75.0 


73.5 


0.8 


82 


25 


242 


IRS13W 


-3.92 


— 2.51 


4.66 


K3III 


-40.0 




-40.0 


1.1 


-74 


30 


244 




4.09 


-2.35 


4.72 


MOIII 


40.9 


39.7 


40.3 


0.8 






398 




1.41 


4.76 


4.97 


MOIII 


— 262.8 




-262.8 


1.1 






416 


IRS7 


-0.26 


5.65 


5.65 


>M7III' 


-108.6 


-109.4 


-109.0 


0.8 


-103 


15 


170 


IRS20 


-0.40 


-5.68 


5.69 


M2III 


-90.7 


-91.1 


-90.9 


0.8 


17 


25 


153 


G577A? 


0.19 


-6.45 


6.46 


M1III 


27.2 


29.8 


28.5 


0.8 


32 


25 


418 




-2.73 


5.92 


6.52 


K3III 


-115.5 




-115.5 


1.1 






189 




4.72 


-4.81 


6.74 


K2III 




26.9 


26.9 


1.1 






308 




6.84 


0.46 


6.85 


K0III 


63.1 




63.1 


1.1 






213 




-5.74 


-3.78 


6.87 


K1III 


41.1 




41.1 


1.1 






144 




-1.20 


-6.96 


7.06 


MOIII 


-23.4 


-21.2 


-22.3 


0.8 






415 




4.45 


5.63 


7.17 


K1III 


68.7 


67.7 


68.2 


0.8 






430 




2.99 


6.57 


7.22 


G8III 




16.0 


16.0 


1.1 






325 


IRS1NB(3) 


7.34 


1.23 


7.44 


M1III 




72.7 


72.7 


1.1 


29 


25 


297 




7.48 


0.06 


7.48 


G9III 




71.1 


71.1 


1.1 






346 


IRS1NE(2) 


7.42 


2.20 


7.74 


M5III 




-68.8 


-68.8 


1.1 


-8 


25 


126 




0.07 


-7.81 


7.81 


MOIII 


-145.3 


-147.0 


-146.2 


0.8 






137 


IRS12N 


-2.81 


-7.50 


8.01 


MOIII 


-60.4 


-60.9 


-60.7 


0.8 


-96 


20 


407 


BHA4E 


-6.08 


5.36 


8.11 


M4III 


-28.2 




-28.2 


1.1 


-77 


30 


459 


G734? 


-3.16 


7.72 


8.34 


K3III 


-28.3 




-28.3 


1.1 


104 


40 


342 


IRSlNE(l) 


8.12 


1.99 


8.36 


M3III 


-67.5 




-67.5 


1.1 


112 


25 


124 


G849? 


3.30 


-8.02 


8.67 


K2III 


-107.3 




-107.3 


1.1 


88 


25 


115 


IRS14N 


1.58 


-8.56 


8.70 


M4III 


-13.9 


-13.3 


-13.6 


0.8 


19 


20 


417 


BHA4W 


-6.89 


5.65 


8.91 


K2III 


88.4 




88.4 


1.1 


100 


30 


467 




-2.62 


8.56 


8.95 


M4III 


— 15.5 


-12.1 


-13.8 


0.8 






151 


IRS9 


6.47 


-6.42 


9.11 


M3III 


-341.4 


-341.3 


-341.3 


0.8 


-300 


25 


96 


IRS14SW 


0.41 


-9.47 


9.48 


M5III 


24.5 


24.6 


24.5 


0.8 


29 


20 


493 


G904 


0.27 


9.51 


9.51 


M2III 


43.4 


44.2 


43.8 


0.8 


23 


30 


101 


IRS12S 


-2.87 


-9.26 


9.70 


M1III 


48.9 


48.7 


48.8 


0.8 


51 


25 


271 


G986 


-9.63 


-1.25 


9.71 


MOIII 


62.7 




62.7 


1.1 


-50 


30 


87 




-1.23 


-9.76 


9.84 


M1III 


-46.9 


-45.1 


-46.0 


0.8 






498 




-0.95 


10.01 


10.06 


K5III 


-5.2 


-1.9 


-3.5 


0.8 






393 


IRS10EE 


8.94 


4.65 


10.08 


MOIII 


-113.1 




-113.1 


1.1 


-55 


30 


290 




-10.19 


-0.28 


10.20 


M2III 


-25.9 




-25.9 


1.1 






94 




-3.80 


-9.60 


10.32 


K4III 


-7.5 




-7.5 


1.1 






174 


G1060? 


8.89 


-5.44 


10.42 


K1III 


25.6 




25.6 


1.1 


-103 


35 


495 


G1044 


-3.89 


9.75 


10.50 


M2III 


199.5 


200.3 


199.9 


0.8 


134 


25 


460 




-8.07 


7.76 


11.20 


MOIII 


-100.6 




-100.6 


1.1 






490 


G1130 


-6.25 


9.35 


11.24 


M3III 


-20.6 




-20.6 


1.1 


7 


30 


338 


G1061? 


11.14 


1.77 


11.28 


M1III 


184.8 




184.8 


1.1 


227 


25 


410 




9.92 


5.48 


11.33 


K2III 


-107.1 




-107.1 


1.1 






171 




-9.90 


-5.56 


11.35 


K3III 


24.5 




24.5 


1.1 






134 




8.57 


-7.51 


11.39 


G9III 


-9.3 




-9.3 


1.1 






437 


G1150 


-9.10 


6.87 


11.40 


M4III 


15.6 




15.6 


1.1 


59 


30 


366 




-11.03 


3.02 


11.43 


K4III 


-32.9 




-32.9 


1.1 






534 




2.02 


11.63 


11.81 


K3III 


-240.6 




-240.6 


1.1 






538 


IRS15N 


0.83 


11.78 


11.81 


M2III 


-12.3 


-10.1 


-11.2 


0.8 


17 


25 


527 


G1152 


-3.42 


11.47 


11.97 


M3III 


-43.5 




-43.5 


1.1 


-68 


30 


68 




-5.35 


-10.71 


11.97 


MOIII 


-64.8 




-64.8 


1.1 






441 


G1174? 


9.70 


7.15 


12.05 


K1III 


-0.1 




-0.1 


1.1 


12 


50 


295 


G1124? 


12.06 


0.00 


12.06 


<G8III 


224.8 




224.8 


1.1 


177 


30 


71 


G1248 


6.45 


-10.23 


12.09 


<G8III 




80.6 


80.6 


1.1 


171 


35 


474 


GUIS? 


8.26 


8.88 


12.13 


K1III 


47.9 




47.9 


1.1 


105 


35 


541 


G1183 


-1.10 


12.09 


12.14 


M2III 


-4.6 




-4.6 


1.1 


45 


20 


404 




-11.16 


5.10 


12.28 


MOIII 


-77.1 




-77.1 


1.1 






464 


G1266 


-9.30 


8.10 


12.33 


K4III 


-113.7 




-113.7 


1.1 


42 


50 


524 


G1198 


-4.98 


11.31 


12.36 


M2III 


-133.1 




-133.1 


1.1 


-88 


50 


491 




-8.19 


9.38 


12.45 


K3III 


38.2 




38.2 


1.1 
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Table 4 — Continued 



Identifier Designation 3 ARA b ADEC b r b Sp. c Vi d V 2 ° V AVE ' CT^ V Gonzcl h <r Gonzol 

Type kms~ 1 kms^ 1 kms^ 1 kms^ 1 kms^ 1 kms -1 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


261 




12.57 


-1.57 


12.67 


G9III 


158.8 




158.8 


1.1 






173 


IRS28 


11.48 


-5.54 


12.74 


M1III 


-47.6 




-47.6 


1.1 


-93 


25 


392 


G1202? 


12.05 


4.51 


12.86 


M3III 


79.9 




79.9 


1.1 


137 


25 


516 


G1187? 


7.05 


10.81 


12.90 


MOIII 


70.4 


71.6 


71.0 


0.8 


37 


30 


528 


G1306 


-6.78 


11.41 


13.27 


MOIII 


11.2 




11.2 


1.1 


-23 


50 


547 


G1311 


5.59 


12.33 


13.54 


K4III 


84.7 


84.4 


84.5 


0.8 


111 


30 


158 


G1390? 


12.18 


-6.32 


13.72 


K4III 


-218.6 




-218.6 


1.1 


-163 


30 


102 


G1438? 


-10.41 


-9.10 


13.83 


M4III 


53.9 




53.9 


1.1 


72 


25 


431 


G1271? 


12.28 


6.77 


14.03 


K1III 


-53.3 




-53.3 


1.1 


114 


25 


506 




-9.84 


10.18 


14.16 


K2III 


-27.7 




-27.7 


1.1 






482 


G1458 


-10.98 


9.12 


14.27 


M4III 


36.8 




36.8 


1.1 


77 


45 


462 


G1334? 


12.26 


7.96 


14.62 


K2III 


-144.4 




-144.4 


1.1 


67 


40 


109 




12.56 


-8.67 


15.26 


MOIII 


-17.9 




-17.9 


1.1 






570 




7.24 


13.60 


15.40 


K2III 


-38.3 




-38.3 


1.1 






539 


IRS11SW 


-10.39 


11.80 


15.72 


M6III 


-91.6 




-91.6 


1.1 


-24 


30 


81 


G1620? 


12.54 


-9.86 


15.95 


G8III 


78.5 




78.5 


1.1 


81 


35 


575 




8.16 


13.75 


15.99 


K3III 


143.3 




143.3 


1.1 






573 


IRS11NB 


-9.31 


13.56 


16.45 


M6III 


-2.7 




-2.7 


1.1 


14 


30 


566 




10.87 


13.26 


17.15 


MOIII 


44.1 




44.1 


1.1 






580 




11.57 


13.86 


18.05 


G9III 


182.8 




182.8 


1.1 







a Designations are taken from Genzel et al. (2000). The GNNNN designations are formed by converting the radius in Genzel ct al. (2000) into a 
string without the decimal point, i.e. G288 is the star with r— 2^88, and G213B is the second star with r— 2''l3. 

b Positions are with respect to Sgr A* and have an error of ±0.1". A negative offset indicates that an object is to the south or the west of 
Sgr A*. 

c Late-type spectral types were determined by comparing spectra to those of stars in the WH atlas (Wallace &z Hinkle 1996). 
d Heliocentric velocities measured from data in first slit-scan. 
°Hclioccntric velocities measured from data in second slit-scan. 
t (V 1 +V 2 )/2. 

s Estimated error in average velocity. 

Heliocentric velocities from Genzel ct al. (2000). 

^RS? has a spectral type of M2I, according to Scllgrcn, Hall, Klcinmann. & Scoville (1987), based on the depth of the CO absorption features, 
the relative lack of water absorption, and the estimated luminosity. Our code only considers CO absorption depth and assigns a very late red giant 
spectral type. 



Note. — Negative velocities indicate a blue-shifted spectrum. 
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Table 5. Velocity Statistics 



Statistic 


GC1 


GC2 


GC1 or GC2 


GC1 and GC2 


(1) 


(2) 


(3) 


(4) 


(5) 


number 


153 


68 






number of cool stars 


78 


29 


85 


22 


average velocity 


-15.37 


9.04 


-10.09 


-3.58 


median velocity 


-12.32 


26.86 


-7.48 


24.52 


standard deviation 


102.18 


100.44 


100.88 


107.65 


average V G ci - V G C2 








-0.80 


< - r V _ GCl-VGC2 








1.53 




1.08 


1.08 




0.77 



Note. — Negative velocities indicate a blue-shifted spectrum. V GC i — V GC 2 is the 
velocity difference for stars that are common in both data sets. The sample of cool 
stars only includes those for which high quality spectra were obtained. 
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Fig. 1. — (a) An image taken with the slit-viewing camera, with synthetic slit apertures 
overlayed to correspond to slit positions used to obtain the first night of data. The slit 
apertures cover 24" in length and (X'72 in width, (b ) The same for the second night of data, 
except that the slit width is 0'.'43. This figure can be found in the full version of the paper 
at http:/ /www-int. stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 2. — Spectral image obtained on June 4, 1999. The image has been processed to 
remove sky and telescope background, dark current, bad pixels, and pixel-to-pixel variations 
in response. North is down. Table 2 gives the approximate wavelength coverage of each 
complete order, with order 33 being near the top of the image and order 38 being near the 
bottom of the image. For reference, note the strong Br-7 (2.166 /xm) emission from ionized 
gas located near the center of the image. The CO bandhead begins in the first complete order 
from the top, about 1/3 of the size of the image from the left. This figure can be found in 
the full version of the paper at http:/ /www- int. stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 3. — SCAM image of the slit-scan region with numerical identifiers (Genzel et al. 
2000) for each star in Table 3. Note the shadow of the slit running approximately 
north-south. This figure can be found in the full version of the paper at http://www- 
int.stsci.edu/ ~figer/papers/nirspec/vel/ms.ps. 
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Fig. 4. — Rectified order from data displayed in Figure 2. The data have been flipped so that 
North is up, and labels identify objects in Table 4. Note that the hot stars (without evidence 
of CO bandhead absorption) are not identified in the table. This figure can be found in the 
full version of the paper at http://www-int.stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 5. — Full spectrum of IRS7 over all echelle orders. This figure can be found in the full 
version of the paper at http://www-int.stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 6. — Spectrum of IRS7 near the CO bandhead feature, as observed (thin), and shifted 
(dotted) by the velocity corresponding to the result from the cross-correlation with respect 
to the spectrum of Arcturus (thick). The spectra are shifted by arbitrary amounts along 
the vertical axis for presentation purposes. This figer can be found in the full version of the 
paper at http:/ /www-int. stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 7. — Spectrum of IRS9 near the CO bandhead feature, as observed (thin), and shifted 
(dotted) by the velocity corresponding to the result from the cross-correlation with respect 
to the spectrum of Arcturus (thick). The spectra are shifted by arbitrary amounts along 
the vertical axis for presentation purposes. This figer can be found in the full version of the 
paper at http:/ /www-int. stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 8. — Spectrum of IRS33E near the CO bandhead feature, as observed (thin), and shifted 
(dotted) by the velocity corresponding to the result from the cross-correlation with respect 
to the spectrum of Arcturus (thick). The spectra are shifted by arbitrary amounts along 
the vertical axis for presentation purposes. This figer can be found in the full version of the 
paper at http:/ /www- int. stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 9. — Spectrum of IRS20 near the CO bandhead feature, as observed (thin), and shifted 
(dotted) by the velocity corresponding to the result from the cross-correlation with respect 
to the spectrum of Arcturus (thick). The spectra are shifted by arbitrary amounts along 
the vertical axis for presentation purposes. This figer can be found in the full version of the 
paper at http:/ /www-int. stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 10. — Spectrum of IRS12N near the CO bandhead feature, as observed (thin), and 
shifted (dotted) by the velocity corresponding to the result from the cross-correlation with 
respect to the spectrum of Arcturus (thick). The spectra are shifted by arbitrary amounts 
along the vertical axis for presentation purposes. This figer can be found in the full version 
of the paper at http://www-int.stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 11. — (a) Relationship between the velocity correction for spectral type versus tem- 
perature for template stars, (b) Relationship between temperature and CO index for 
template stars. This figer can be found in the full version of the paper at http://www- 
int .stsci.edu / ~figer / papers / nirspec / vel /ms.ps. 
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Fig. 12. — (a) Histogram of velocities for cool stars in the central parsec from Table 4. (b) 
Same, using data from Genzel et al. (2000). This figer can be found in the full version of the 
paper at http:/ /www- int. stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 13. — Plot showing radial velocities for cool stars in the central parsec from Ta- 
ble 4. Open (filled) circles correspond to stars with motion away from (toward) the 
Sun. The area of each circle scales with velocity. The origin corresponds to the posi- 
tion of Sgr A*. This figer can be found in the full version of the paper at http://www- 
int.stsci.edu/ ~figer /pap er s/nir sp ec / vel/ms . ps . 
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Fig. 14. — Comparison of heliocentric velocities for the 49 stars in common between our 
sample and that in Genzel et al. (2000). Errors for the Genzel et al. data are shown with 
error bars. Errors for the data in Table 4 are smaller than the symbol size. A line with 
a slope of unity is overplotted. This figer can be found in the full version of the paper at 
http : // www-int . stsci.edu / ~figer / papers / nirspec / vel/ms. ps . 
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Fig. 15. — (a) Radial velocity of stars in the sample as a function of galactic longi- 
tude and (b) galactic latitude. The horizontal axes are projected distances from Sgr A* 
in parsecs. This figer can be found in the full version of the paper at http://www- 
int .stsci.edu / ~figer / papers / nirspec / vel /ms.ps. 
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Fig. 16. — (a) The thick solid line represents the normalized asymmetry statistic, A/N, 
versus position angle, from 0° to 180° East of North, as calculated using equation I. A 
is the number of stars with negative velocities minus those with positive velocities on one 
side of the axis minus the same quantity for stars on the other side of the axis, and N is 
the total number of stars in the data set. The dotted and dashed lines show the number 
of stars with negative velocities minus those with positive velocities on each of the two 
sides of the axis, divided by N. The plots with the highest values show the total number 
of stars on each side of the axis, (b) Cumulative percentage of systems versus maximum 
normalized A for 1000 simulated systems and all position angles. Each system has 85 stars 
with uniform random location within a 1 pc square area centered on Sgr A* and a gaussian 
velocity distribution having mean of zero and standard deviation of 102 kms~ l . Note that 
20% of systems have A >0.2. This figer can be found in the full version of the paper at 
http : // www-int . stsci.edu / ~figer / papers / nirspec / vel/ms. ps. 



Fig. 17. — Same as Figure 16, except data from Genzel et al. (2000) have been used. The 
Monte Carlo simulation uses systems with 198 points, the number of late-type stars in the 
Genzel et al. (2000) data set. This figer can be found in the full version of the paper at 
http : // www-int . stsci.edu / ~figer / papers / nirspec / vel/ms. ps. 
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Fig. 18. — (a) Velocity dispersion in 5" wide annuli as a function of radius, (b) Same, using 
data from Genzel et al. (2000). This figer can be found in the full version of the paper at 
http : // www-int . stsci.edu / ~figer / papers / nirspec / vel/ms. ps. 
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Fig. 19. — (a) Enclosed mass according to equation 2 (Bahcall & Tremaine 1981) using 
velocities adjusted to the Local Standard of Rest by adding +10 kms -1 , as a function of 
radius, (b) Same, using data from Genzel et al. (2000), after subtracting the median velocity 
of the sample in order to adjust the data to the LSR. IRS9 has been excluded from the 
analysis for both data sets because of its anamolously high velocity. This figer can be found 
in the full version of the paper at http:/ /www-int. stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 20. — (a) Surface number density in our data, (b) Surface number density in Genzel et 
al. (2000) data. Neither plot has been adjusted for incompleteness. 



Fig. 21. — (a) Observed and model velocity dispersion in the central parsec. The model 
assumes circular velocities, a central dark mass of 3.45(10 6 ) M , particles with non-zero 
mass, and a stellar density distribution according to equation 3 with r =0.4 pc and n=— 0.25 
(solid) and n=— 2 (dashed), (b) Same, using data in Genzel et al. (2000). 



Fig. 22. — (a) x 2 for model compared to observations of the velocity dispersion as a function 
of power-law index for stellar density according to equation 3 with r =0.4 pc. (b) Same, 
using data in Genzel et al. (2000). 



Fig. 23. — (a) Locations of 10,000 particles in simulation for 1/r 2 distribution. 
(b) Same, assuming a distribution according to equation 3 with r =0.4 pc and 
n=— 0.25. This figure can be found in the full version of the paper at http://www- 
int.stsci.edu/ ~figer /pap ers/nir sp ec / vel/ms . ps . 
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Fig. 24. — (a) Projected number density (left axis) and velocity dispersion (right axis) of 
model system, assuming a 1/r 2 distribution at all radii, (b) Same, assuming a distribution 
according to equation 3 with r =0.4 pc and n=— 0.25. This figure can be found in the full 
version of the paper at http://www-int.stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 25. — (a) Enclosed mass in the Monte Carlo simulation and that inferred by using 
equation 2 (Bahcall & Tremaine 1981), assuming circular orbits, as a function of radius, 
and a system having a 1/r 2 distribution, (b) Same, assuming a distribution according to 
equation 3 with r =0.4 pc and n=— 0.25. This figer can be found in the full version of the 
paper at http:/ /www- int. stsci.edu/~figer/papers/nirspec/vel/ms.ps. 
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Fig. 26. — Line-of-sight acceleration versus projected distance from Sgr A* for a star 
on a circular orbit. We have assumed a central point mass of 3.45(10 6 ) M , an or- 
bital phase of 7r/4, and a range of inclinations between the orbital axis and the line 
of sight. This figure can be found in the full version of the paper at http://www- 
int .stsci.edu / ~figer / papers / nirspec / vel /ms.ps. 



